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Abstract 

The level-spacing distribution of consecutive energy eigenvalues is calculated numerically 
at the metal insulator transition for 3d systems with different cuboid shapes. It is found 
that the scale independent critical P c (s) changes as a function of the aspect ratio of the 
samples while the critical disorder W c /V =16.4 remains the same. We use our data to test 
whether an expression for the small-s behaviour of the level statistics proposed by Kravtsov 
and Mirlin for the metallic regime is applicable also at the critical point. For this reason, a 
shape dependent dimensionless critical conductance g c has been extracted from the small-s 
behaviour of the critical level statistics. Our result for a cubic sample, g c — 0.112±0.005, is in 
good agreement with a value obtained previously from calculations using the Kubo-formula. 



1 Introduction 

Energy eigenvalue correlations provide general tools for the statistical description of disordered 
materials. Among them, the nearest neighbour level spacing distribution P(s) represents one of 
the simplest statistics. Here, s = \Ei + \ — Ei\/A is the energy difference of consecutive eigenvalues 
Ei divided by the mean level spacing A. Nevertheless, knowing P(s), it is possible to tell whether 
a given system exhibits metallic or insulating behaviour at temperature T = K. 

Considerable attention has recently focused on the special case of critical statistics which 
were found at the metal insulator transition in disordered 3d [||, |2|, || ||, ^, ^, ||, |j| and 2d 



p0| ) [T2j , 13] systems. In all cases, the numerically obtained critical level statistics were noticed 
to be independent of the system size L. This is different in the metallic and insulating regimes 
where a size dependence has been observed. In the limit L — > oo, the nearest neighbour spacing 
distribution follows the Poissonian decay, P(s) = exp(— s), for disorder strength W larger than 
the critical disorder W c . The metallic side, W/W c < 1, is well described by random matrix 



theory (RMT) JbJ, |15|, |16|, |17| . As in the metallic phase, the critical P c (s) depends on the 
symmetry of the Hamiltonian describing the system under consideration. In addition, and in 
contrast to the universal RMT description for the metallic phase, the critical P c (s) depends also 
on the spatial dimension. The respective forms of P c (s) are known only from numerical studies. 

In Ref. [J/J using the Anderson model, the application of an Aharonov-Bohm-flux was shown 
not only to change the universality class of the level statistics due to breaking of time reversal 
symmetry. It has also been demonstrated that, depending on the strength of the flux, there exists 
a set of scale independent critical spacing distributions. These are associated with a continuous 
crossover from orthogonal to unitary symmetry |7|, [l8| . It is well known that the corresponding 
Hamilton matrix can be transformed such that the AB-flux is completely absorbed into the 
boundary conditions. So the question arises, whether or not P c (s) does in general depend on the 
boundary conditions applied. Recently, it has been shown |]l9| that P c (s) is in fact sensitive to a 
change of the periodic boundary conditions, which were always applied previously, to Dirichlet 
boundary conditions in one, two and three directions. 
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Figure 1: The critical level spacing distribution P c (s,q) versus the level spacing s for aspect 
ratios q = 1 Lq = L z = 30 a), q = 1/7 (o, Lq = 56 a, L z = 8 a), and q = 20 (•, Lq = 
10a,L z = 200a). The solid line is the Poissonian decay exp(— s). The insert displays the 
semi-logarithmic plot of the large-s behaviour. 

In this paper we show that P c (s) is sensitive to the shape of the sample, too. We present 
results of a numerical investigation on various finite size cuboids with square base Lq and height 
L z , characterized by the aspect ratio q = L z /Lq. Assuming that a formula for the small-s 
behaviour of P(s), proposed by Kravtsov and Mirlin [^0| for the metallic regime, may approxi- 
mately be used at the metal-insulator transition, a dimensionless critical conductance g c can be 
extracted from the small-s behaviour of the calculated P c (s). We find that g c also depends on 
the sample shape and that the value obtained for the cubic system is in close agreement with a 
result obtained by a different method. 



We investigate the conventional 3d Anderson model which is described by the Hamiltonian 



The creation, (4, and annihilation, c m , operators act on the states of non-interacting electrons 
at the sites {to, n} of a 3d simple cubic lattice. The disorder energies e m are chosen to be a 
set of independent random numbers distributed in the interval [—W/2,W/2] with probability 
1/W, where W denotes the disorder strength measured in units of V . We consider only transfer 
between nearest neighbour sites. Periodic boundary conditions are applied in all directions. The 
size of the different systems investigated are described by a quadratic base Lq and a height L z 
ranging from 8 to 200 lattice constants a. The shape of the system is conveniently denoted by 
the dimensionless aspect ratio q = L z /Lq. For all shapes, a size independent metal insulator 
transition is found at a critical disorder of W c /V = 16.4. 

The eigenvalues used in our study have been taken from an interval [—0.5 V, 0.5 V] around 
the band centre. They are obtained by direct diagonalization using a Lanczos algorithm. For 
each aspect ratio, several system sizes are computed with a large number of different realizations 
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Figure 2: The second moment of the critical level spacings, I c (q), versus the aspect ratio q = 
L z j Lo-The corresponding system sizes are in the range from (Lq/o, = 80, L z /a = 8) to (Lo/a = 
W,L z /a = 200). 

of the disorder potentials so that for each particular case the number of eigenvalues exceeded 
10 5 . 



3 Results and Discussion 

In Fig. H| the size independent critical level spacing distribution P c (s, q) is shown for different 
aspect ratios: a flat quasi 2d sample (q = 1/7), a cube (q = 1), and a long narrow bar (q = 20). 
Deviations of the aspect ratio from q = 1, result in critical spacing distributions that seem to 
approach the Poissonian decay, P(s) = exp(— s), as q — > (quasi 2d sample) or q — ► oo (quasi 
Id sample). 

To give a more quantitative description, we have calculated the second moment I c (q) of the 
critical level spacings 

poo 

I c (q) = 1/2 < s 2 >= 1/2 / s 2 P c (s, q) ds (2) 

Jo 

which is shown in Fig. [2] as a function of the aspect ratio q = L z /Lq. The cube (q = 1) represents 
the minimum value of the second moment, I c (l) = 0.705 ± 0.005, whereas I c (q) — * 1 for aspect 
ratios that strongly deviate from the cube. Our value of I c {q = 1) is in accord with the result 
of Ref . H . From the symmetric appearance of the semi-logarithmic plot in Fig. ^ it becomes 
apparent that I c (q) — * Ic0-/q) when crossing q = 1. 

At present, there exists no complete analytical theory which describes the critical energy 
level correlations. To compare our results, we therefore have to put up with a formula proposed 



recently |2C] for the metallic regime. The two-level correlation function is defined as 
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R(s) 



<p(E)p{E + uj) >, 



(3) 



< p > z 

where p{E) is the density of states at energy E, s = uj/A, and < ... > denotes averaging over 
realizations of the disorder potential. In the metallic regime and for small s, R(s) becomes |^] 
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Figure 3: The shape dependent critical conductance g c (q) as a function of the aspect ratio q. 



This expression contains corrections to the usual Wigner-Dyson form in the region uj ~ A << E c , 
where E c is the Thouless energy. Eq. (|j) depends on the dimensionless conductance g >> 1, 
the spatial dimensions, the boundary conditions and also on the shape of the system via the 
diffusive modes [p0]. For cuboids L\ x L z and periodic boundary conditions one gets 
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At the critical point, where g = g c < 1, the result of the expansion made in Ref. [20] is expected 
to give only qualitatively correct behaviour. In order to test this expectation, we tentatively 
assume here that Eq. (Q) actually holds at the metal insulator transition. Therewith, we are able 
to extract a shape dependence of the size independent critical conductance g c from the small-s 
part of the calculated P c (s,q), because for small s, R(s) = Y2k s) + S(s) ~ P(0,s), where 
P(k, s) is the probability density to find exactly k eigenvalues within the interval s. 

The critical conductance as a function of the sample shape is shown in Fig. [|. Here, the 
values for q = 0.1 and q = 20 have been omitted, because it was not possible to accurately 
extract the small-s gradient from the data. Again, the already noticed symmetry in replacing 
q — * 1/q when crossing q = 1 becomes evident. 

Our result for the cubic system, g c (l) = 0.112 ± 0.005, agrees very well with the value 
g c = 0.10 ± 0.01 obtained from the numerical calculation of the Kubo-conductivity [21|. Since 
the two methods for determining the critical conductance are completely different, the perfect 
agreement of the numerical values is remarkable. If this accordance is not mere accidental, it 



has to be understood why the result of the expansion in Ref. [ 20 ] holds also quantitatively at the 
critical point. For Dirichlet boundary conditions, we obtain a numerical value g c = 0.079 ±0.002 
for the cubic system. 

In mesoscopic systems, the decay of g c {q) accompanied with deviations from the cubic shape 
may explain some of the different values (0.03-0.2) p^] reported for measurements on various 
samples. However, the question remains to be answered, whether or not the expression for 
the small-s behaviour proposed by Kravtsov and Mirlin |^] for the metallic regime is really 
applicable at the metal insulator transition. 
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In conclusion, a dependence of the scale independent critical level statistics on the shape of 
the samples has been detected at the Anderson transition in 3d systems. Using a formula which 
has been proposed for the metallic regime [^(J to fit our data obtained at the critical point, we 
find a shape dependence of the scale independent conductance g c . Provided that this method 
of extracting the conductance is justified also at the critical point, the various experimentally 
obtained critical conductances reported in the literature can possibly be attributed to different 
shapes of the samples investigated. 
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